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Abstract: Cosmic superstrings are produced towards the end of the brane inflation, a 
scenario realized in modern superstring theory. If the string tension is low enough, loops 
tend to be relatively long-lived. The resultant string network is expected to contain 
many loops which are smaller than typical Galactic scales. Cosmic expansion damps the 
center of mass motion of the loops which then cluster like cold dark matter, eventually 
decaying by emission of gravitational radiation. 

Loops will lens stars within the galaxy and local group. We explore microlensing of stars 
as a tool to detect and to characterize some of the fundamental string and string network 
properties, including the dimensionless string tension Gfi/c^ and the density of string 
loops within the Galaxy. 

As G/i the intrinsic microlensing rate diverges as 1/ ^/GJ^ but experimental detection 
will be limited by shortness of the lensing timescale and/or smallness of the bending angle 
which each vary oc Gfi. We flnd that detection is feasible for a range of tensions. As 
an illustration, the planned optical astrometric survey mission, Gaia, should be able to 
detect numerous microlensing events for string networks with tensions 10~^°^G';U< 10~^. 
A null result for optical microlensing implies Gfi ^ 10~^°. 

If lensing of a given source is observed it will repeat because the internal motions of the 
loop are relativistic but the center of mass motion may be much smaller, of order the 
halo velocity Vh- This distinctive hallmark, c/vh ~ 1000 repetitions, suggests a useful 
method for confirmation of a potential lensing detection. 

We argue that the estimate of the Galactic lensing frequency is likely to rise with more 
realistic descriptions of the superstring network while the effect of the inclusion of loop- 
loop interactions within the Galactic halo is not yet clear. 
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1. Introduction 

Inflation probably explains the origin of our universe A concrete realization of this 
scenario in superstring theory is brane inflation [0, ^, ^. The simplest version involves 
an interacting D3-D3-brane system in which the D3-brane moves towards the D3-brane 
sitting at the bottom of a warped throat, with the position of the mobile D3-brane 
corresponding to the inflaton. Inflation ends when the D3-brane and the D3-brane collide 
and annihilate, initiating the hot big bang. Cosmic strings are copiously produced during 
the epoch of annihilation ||5|, Q . These strings are simply superstrings (F- and D-strings 
and their bound states 0, ^, |^) stretched to cosmological sizes. Finding evidence for the 
existence of these strings would go a long way towards conflrming certain fundamental 



aspects of superstring theory |I0|, |Tl 



The single most important property of a cosmic string is its tension fi, or, in dimen- 
sionless terms, the characteristic gravitational potential Gfi/c^ [jl2[ (hereafter, we will 
take c = 1). An initial estimate of the range of tension produced towards the end of 
inflation gave 10^^^ ^Gfi^ 10^^ 0. A recent analysis of the reheating that accompanies 
brane inflation, including multi-brane multi-throat scenarios, suggests that Gfi < 10^^^ 
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is quite possible [|^. The properties of the cosmic superstrings produced at the end of 



inflation are compatible with all observations today but are beginning to be constrained 



by observational data |T^, |T5|, |1^]. The actual bound is sensitive to the details of the 
string tension spectrum and the probability of string-string interactions but Gfi ^ 10~^ 
gives a rough idea of current observational constraints. 

For small Gfi, it will be very challenging to detect cosmic string signatures in lens- 
ing, gravitational wave bursts, pulsar timing and B mode polarization in the cosmic mi- 
crowave background radiation by methods that have been investigated hitherto. How- 
ever, even low Gfi strings produce double images of point-like objects. Here, we will 
investigate the utility of microlensing of stars to learn about the string content of the 
universe. The propensity to microlens depends, of course, upon the number of lenses 
and their cross section for bending light. Our current understanding of string network 
evolution implies that the rate of lensing increases as Gfi decreases. This happy circum- 
stance may eventually lead to practical experiments to detect superstrings. Microlensing 
of stars is much more promising than microlensing of quasars [p!7| . 

Although it is well accepted that cosmic strings evolve to a scaling network |I8| , 
string network evolution is by no means well understood. Even if one had a complete 
knowledge of the intrinsic string tensions in string theory, the cosmological evolution of 
the string network remains a challenging problem [|ll] . Recent analyses strongly suggest 
that cosmic superstrings evolve dynamically so as to produce a scaling solution in which 
there exists a stable relative distribution of strings with different quantum numbers 



IS, 20 1 . This is very much like the behavior of cosmic strings generated by abelian 



Higgs or Nambu-Goto type models |18, 12|. The superstring solutions found in current 
string theory, however, form a much larger class which is only now being explored. 

A scaling solution implies the fractions of the critical energy density in long strings 
and loops are constant. In fact, the determination of Qiong and Qioop has been an out- 



standing question for some time |T2[. Numerical simulations for cosmic strings may be 
fitted by 

Qs = ^long + ^loop = i'Gfl + X^Gji (1.1) 

for constants and x- A small isolated loop of length / will decay in a characteristic time 
T = l/{rGfi) and F ~ 50 for Nambu-Goto strings. A common assumption made in many 
simulations is that small loops are produced by the intersection of long strings (including 
self intersection). For Ht << 1 where H is the Hubble constant the loops decay quickly 
via gravitational radiation and Qg ~ ^long ~ ipGfi; the loop contribution is negligible. 
Numerical simulations show ~ F. Presumably, this is a good approximation for G/j, 
close to the present day observational bound. 

Several new considerations have emerged from the latest simulations and from the- 
oretical work. 

• Loops dominate: For small Gn newly formed loops do not immediately dissipate. 
The latest simulations ETI, ESl E3] imply a loop distribution which, if truncated by 
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gravitational wave damping, yields a qualitatively new effect: ~ ^loop » ^long- 



Small scales on long strings: Analytic work [T^ shows that the stretching 
of long strings is insufficient to prevent the development of small scale structure. 
Collisions of long strings may produce large numbers of loops of all sizes. Estimates 
indicate ~ 10% of the long string length ends up in loops with scale comparable 
to the horizon. Loops with Ht > 1 contribute to Qs- 

Backreaction of small scales: When Qioop increases, small loops can recombine 
to form bigger loops and, likewise, isolated loops may be incorporated into long 
strings. The determination of the scaling solution when energy can flow up and 
down a hierarchy of scales is an outstanding problem. 



• Damping of loops vs long strings: Analytic work p5| , ^ indicates that the rate of 
damping of modes on long strings differs from TGfi, the characteristic loss rate for 
an isolated loop. Under such circumstances the network- and time-averaged gravi- 
tational energy loss rate per mode is coupled to the details of the intercommutation 
process. 

These factors motivate the development of a very approximate description of the 
string network when the tension is low. Specifically, we extend the analysis of Olum 



et al. in Ref.[^ to describe the loop population for small Gjj. We also parameterize 
the effective decay rate of loops in the network by TuGfi where, in general, we expect 
F/j < r ~ 50. 

This minimal quantitative framework is used to begin an investigation of microlens- 
ing by loops. Lensing by cosmic string loops at cosmological distances was considered 
originally in |32|. The analyses focused on Gfi ~ 10^^ and the possibility that the 
lensing might produce multiple images of QSOs. The situation of interest here is quite 
different though the physical process, gravitational lensing, remains the same as was 
originally discussed in these papers. We shall first give a brief sketch of the latest sta- 
tus of the cosmic string network. We shall outline the form of the loop distribution 
which is presented in full detail in the Appendix. We shall then use the result to esti- 
mate the event rates and the detectability of microlensing for the cosmic strings. While 
many properties of superstrings and astrophysics remain unaccounted for (specifically, 
the intercommutation probability and the collisional interactions within the Galaxy) this 
approach highlights the role of lowering Gfi and sets a benchmark that we will refine as 
the input physics is better understood. 

Our analysis shows that the intrinsic rate of microlensing increases as tension de- 
creases. At the same time the detection efficiency for a microlensing event diminishes 
because both the duration of the event and also the size of the deficit angle diminish.^ 



^The deficit angle due to tlie string must be bigger tlian the angle subtended by a star. For typical 
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To illustrate the potential utility for detecting cosmic string microlensing of stars, 
we consider the upcoming European Gaia mission |2^. Gaia is primarily an astrometry 



satellite designed to measure proper motions and parallaxes of about 10^ stars in and 
near the Galaxy over a 5 year period. The satellite records the flux for each star and it 
is this feature which is of greatest interest to us. On average each star will be observed 
80 times; individual observations take about ~ 3 seconds. The expected accuracy of the 
photometry from a single observation is easily sufficient to see a factor of 2 change in 
flux from one observation to the next. Using a crude model of detection efficiency, we 
find that an instrument like Gaia would be able to detect microlensing events for cosmic 
strings in the Galaxy: a few events for ~ 10"^'^ to a few dozens for Gfi ~ 10^®. This 
estimate invokes several conservative choices including the value of F/?; the event rate 
increases for < F. 

A key reason for this detectable rate is the over-density of loops in the Galaxy. 
Even small loops can survive a Hubble time when the tension is low. A generally new 
scenario for cosmic string evolution now unfolds: the velocity of the center of mass of 
a loop decreases on account of cosmic drag; loop-loop interactions freeze out. Loops 
subsequently behave like cold dark matter, slowly radiating gravitational waves. Most 
of the loops of interest for microlensing were born during the radiation era and fall into 
gravitational perturbations that begin to grow after equipartition. In particular, the 
loops track the perturbations of the cold dark matter. This clumping enhances the local 
loop density by at least a factor of ~ 10^. 

From the string theory point of view, the true lensing rate may be enhanced above 
our estimate because the intercommutation probability of superstrings may be as small 



as P ~ 10 |T9|. Also, cosmic superstrings come in a variety of tensions and charges 



so that a number of species are present in the network EO]. These effects will tend to 



increase the energy density throughout the universe in the superstring network in eq. 



LTD roughly like fi^ n^lg/P, where n is the effective number of types, n ~ 5. The 
spectrum of cosmic superstrings yields a discrete set of tensions that can easily vary by 
an order of magnitude. 

From an astrophysical point of view, the true lensing rate may be enhanced by 
dissipative loop-halo interactions that boost the galactic loop density and it may be 
diminished by the collisional interaction of loops once they begin to clump in the Galaxy. 
Intercommutation might chop up small loops into smaller loops and shorten their lifetime 
to gravitational wave emission. This, in turn, may diminish the part of the spectrum 
responsible for most of the lensing while increasing the locally generated gravitational 



sources (a solar mass main sequence star within the halo) this is roughly G/i > 10^^"'. Likewise, the 
duration Stiens of the shortest detectable optical microlensing event is limited by, among other things, 
the ability to detect a factor of 2 change in the number of the source photons. If we assume a typical 
source (a solar mass star at distance 10 kpc), a broadband instrument detecting photons near the peak 
of the spectrum (~ 1 eV), a typical instrument (a 10% efhcient meter-class telescope), then we can 
estimate that a 5a detection requires Stiens ^ 40 ms. This corresponds to roughly Gfi > lO"^"*. 
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wave background. We have not included either the dissipation or the interaction in our 
microlensing estimates but will return to these issues in the future. 



2. String Network 

We shall start by considering cosmological cosmic strings like Nambu-Goto strings or 
vortices in an abelian Higgs model. One readily identifies two components, long horizon- 
crossing segments and sub-horizon closed loops. Only straight strings are static; all 
others are dynamic with relativistic motion. Even isolated loops oscillate in a highly 
complicated fashion. A network of long strings and loops can change topology by the 
process of intercommutation, the breaking and rejoining experienced when the motions 
of two segments of string cause them to coincide in 3D space. This process can fragment 
as well as rejoin the basic elements of the string network. Self intersections of a long 
string cuts off new loops, self intersections of a loop transforms it to multiple loops, 
etc. Conversely, inverse processes allow loops to reconnect to long strings and loops to 
reconnect with other loops. In the cosmological context the processes of intercommuta- 
tion, damping (e.g. gravitational wave emission) and cosmological expansion govern the 
string network evolution. 

When a network exhibits scaling behavior the energy density (either long or loop) 
is a fixed fraction of the critical energy density eq. (|1 . 1|) . When very small loops are 
formed and decay promptly the critical density 

Qs ^ ^long ^ i^Gfl (2.1) 

and simulations show ip ^ T where F ~ 50 for abelian strings. The scaling solution 
depends upon F in an indirect way. Small scale structure on the long strings is damped 
by gravitational wave emission with rate cx FG/i if one assumes modes damp like those 
of an isolated loop. An intersection of two long string segments is more dissipative if it 
converts a greater overall length to small loops. In fact, more loops are produced if the 
segment possesses more small scale structure. The dissipation per collision of long string 
segments, therefore, scales oc l/FG/i. The rate for one segment of a horizon-crossing 
string to encounter another oc fi/ong- Therefore, the total dissipation rate per long string 
segment oc Qiong/^Gfi. A fixed rate of dissipation per string is essential to achieve a 
scaling solution in the first place, so that fliong oc TGfi. 

The specific scenario above depends on the Ht << 1 where H is the Hubble constant, 
r = I /TGfi is the damping time and / is the loop size. In a scaling solution a typical loop 
formed at time t has size / ~ at where a is constant so a /TGfi << 1 is required. The 
value of a is poorly known; a < 10~^^ is sometimes invoked but recent numerical 
simulations ||23|, |28| and analytic studies suggest that loops with a range of sizes 



10~^ ^ a ^ 0.25 are created. Lower tension implies longer damping time so that the 
prompt loop decay which is characteristic of the scenario above is now more difficult to 
achieve. 
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To build a model capable of describing the loop distribution for small G/i we follow 
the line of reasoning of Olum et al. [0]. In the radiation or matter-dominated eras the 
number of loops produced per unit loop length per unit volume per unit time has the 
form 

^^^rVW w,th.^/A (2.2) 

for some function f{x) for a scaling solution.^ Numerical simulation in Ref . [p^, |28|1 
suggest a power law distribution of loops 

f{x) = Ax-^ for X < a (2.3) 

where at is the largest loop scale and /? < 2. The constants A, a and (3 may be extracted 
from time-dependent simulations for cosmic strings which appear to have entered a 
scaling regime. 

Assume that in each infinitesimal time interval {t, t + dt) the network produces the 
loops described with / < at according to dN/ dldtdV and these are subsequently diluted 
by cosmological expansion without further intercommutation. This is an instantaneous 
fragmentation description adjusted to agree with the results of time-dependent simula- 
tions that model the complete process. A non-interacting loop shrinks by gravitational 
wave emission until it disappears in time 1/{TrG^). Note the introduction of in the 
dimensionless decay rate in place of F ~ 50. Here, F/j < F implies that loops live longer 
because of complex unmodeled network effects. For example, in the conventional view, 
a loop of size / emits energy at half the total rate of two loops of size 1/2. In addition, 
[p^ , p7| find that the damping rate of long strings is oc (Gfi)'' with k > 1. For the simple 
energy loss recipes employed it makes a great deal of difference whether a length is part 
of a small loop, a large loop or a long, horizon-crossing string. We do not attempt to 
disentangle these effects but simply introduce F^j as an effective damping rate. 

The superstrings differ from previously studied cosmic strings in fundamental char- 
acteristics like Gfi and P and possibly in the resultant loop distribution described by A, 
a and p. It is tricky to extrapolate to describe domains not previously simulated. Our 
calculation of the distribution function for loops is given in detail in Appendix 1. The 
general idea is as follows: When the network scales the long strings are chopped to loops 
whose total length is a constant fraction of the horizon. Energy conservation implies 
that / fxdx is fixed. If the typical loop size a were to be modified (for example, because 
the intercommutation probability is decreased), then this energy conservation argument 
gives A (X Using this reasoning we can, in principle, consider various descriptions 

for /, cosmological expansion dynamics, F/j and Gfi. The minimal assumption is that 
(1 — {v'^))/'^1 is the same for undamped Nambu-Goto strings in Minkowski space as for 
damped superstring networks in FRW cosmology; here, the characteristic inter-string 



^This form does not apply during the recent, A-dominated phase since the horizon is not cx t. The 
foops of direct interest to us are from earUer epochs. 
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distance d{t) = 7st and (f ^) is the square of the string velocity averaged over the length 
of long strings. For our applications, however, we also assume the simulation-derived 
quantities like f3 are fixed and concentrate on changes to the network that arise from 
varying G/i, Tr and the cosmological dynamics. 

At any time the distribution of non-interacting loops in a given volume is the inte- 
grated production rate over the history of the universe 

subject to upper and lower cutoffs. The upper cutoff corresponds to the time when 
gravitational wave damping removes the loop; the lower cutoff is the earliest epoch when 
the loop might be produced. If the expansion scale factor varies a oc then the integrand 
oc jt3n+/3-5 Yoi both radiation and matter dominated eras, the loop distribution at length 
/ is typically dominated by production at early times t = l/a (this requires /? < 7/2 or 
3, respectively) so that 

dN /^"-^ 

dldV ^ t^^a^^-^ ^ ■ ^ 

as previously demonstrated by Vanchurin et al. p8l. This distribution may then be 



compared to the cosmological simulation results ||25|, None of the simulations in- 
corporate gravitational damping so that the cutoff at small loop size in the simulations 
reflects either initial conditions, finite resolution and/or finite simulation times. Nonethe- 
less, the results at various epochs illustrates the buildup of a powerlaw distribution of 
loops at intermediate scales with the predicted slope. In the radiation era, for example, 
eq. i WM ) gives dN/dldV oc while the simulations of Ringeval et al. |^ imply 
dN/dldV oc 

The energy density, lensing probability and lensing rate of loops all involve essentially 
the same moment of the loop distribution, J Idl-^^. The cutoff / = TjiCfxt is the key 
parameter which varies with the tension when we evaluate those quantities. For loops 
originally generated in the radiation era 



fx I a 



1/2 

Since pcr oc we infer 

Sloops oc (2.7) 




This illustrates the combination of parameters that determines fi^oop; if ^long oc G/i this 
square-root behavior suggests that VLioop > ^long for small G/i . 

The numerical simulations which include both loops and long strings allow a quan- 
titative check. Assume that the slope and amplitude of the simulation-derived loop 
distribution is extended to the gravitational cutoff and that the long string density is 
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fixed. The ratio is 



_ J 4.4xlO-y^ Olumetal. 
^loop ) 7.5 X 10-^ (tP^)°'^ Ringeval et al. H 



(2.8) 



which shows the dominance of the loops. The evaluation of the lensing probability and 
rate involves the same integrals which depend on the gravitational wave damping cutoff 
in the same way. When strings have low tension loops, not horizon crossing strings, are 
the favored lensing candidates. 

The appendix treats general powerlaw forms for / with upper and lower cutoffs {au, 
«!,) and various slopes (/?) and, most significantly, various Gfi and radiation loss rates Fr. 
It gives approximate expressions valid for loops from both radiation and matter epochs. 
It compares the simple analytic approximations to more realistic A-CDM cosmologies. 
It provides a range of numerical solutions illustrating how the loop distribution varies 
with fi, Tn, P, and au- 

As mentioned earlier, cosmic superstrings have different properties from cosmic 
strings. The intercommutation probability of vortices is known to be around unity, 
P ^ 1, while that of superstrings is rather complicated, but P ^ g1 |T^, where the 
string coupling g^ ~ 1/10. Also, cosmic superstrings come in a variety of tensions and 
charges so that a number of species are present in the network These effects likely 



increase the energy density in the superstring network compared to its cosmic strings 
counterpart roughly 

n 

a ^ p^s (2.9) 

where n is the effective number of types, n ~ 5. For very small P, it has been argued 
that 1/P 1/P^/^ I^O]. The implication is that the number density of small cosmic 
superstring loops will be boosted with respect to cosmic strings. 



3. Microlensing 

Microlensing refers to brightness variations of a background source caused by a changing 
gravitational field somewhere along the line of sight. The field may be generated by a 
dark point-like object and the astrophysically anticipated candidates include dim stars, 
white dwarfs, neutron stars, and black holes. The paths of photons emitted by the source 
are bent so that the angular area of the source visible to the observer changes. In short, 
the flux from the source changes. The key point is that one does not have to resolve the 
source or the lens to observe the change. Lensing distorts a single source image or, if the 
impact parameter of the photon is within the Einstein radius, lensing creates multiple 
images. The area within the Einstein radius is quite small so its much more common 
for magnification of a single image to occur than the creation of multiple images. If 
source, lens and observer have constant velocities, the time- dependent magnification has 
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an a "priori known functional form. Previous surveys (MACHO, OGLE, EROS, etc.) 
have searched for and identified numerous events with the expected time-dependent 
achromatic form. 

Now, instead of the normal astrophysical lensing candidates, consider a stationary 
straight infinite string oriented perpendicular to the observer's line of sight with respect 
to a background source. Let two photons from the source travel toward the string. 
The photons do not suffer any relative defiection during the fiy-by as long as they pass 
around the string in the same sense. Images formed from photons are undistorted. This 
contrasts with the shear and distortion produced by a point mass. 

However, there is a small angular region just like the Einstein radius that yields 
multiple images as long as the source itself is small in angular size compared to the 
Einstein radius. In essence, some of the photons pass around the the string in a clockwise 
sense and others do so in a counterclockwise sense. Two paths from the source to the 
observer yield two undistorted images. 

Sources that lie behind the string and within the characteristic angle ~ SttG// will 
appear as double images. Unresolved, lensed sources will appear to fiuctuate in brightness 
by a factor of 2 as the angular region associated with the string passes across the observer- 
source line of sight. The characteristic Einstein angle is 



The characteristic angular size of a stellar source at distance R is ©0 = Rq/R. The 
relative size is 



which shows when the stellar source will generally be well described as a point source. 
The relativistically moving and oscillating string will create brightness fiuctuations in 
the background star that can be searched for in a microlensing experiment. 

The actual situation is somewhat more complicated. For a loop, as opposed to a 
straight string, one expects lensing like that of a point mass for photons with impact 
parameter large compared to the size of the loop and lensing like that of a straight string 
for paths that pass close to a segment of the string. We will eschew the complications 
associated with small scale structure on the string and concentrate on photons that pass 
close to a smooth segment of the loop. 

The characteristic scale of the smallest loops today is 




(3.1) 




(3.2) 



Ig = TRG^ttoday = 41pC ( ^^pr j I 13" 




(3.3) 



and the characteristic mass scale associated with such a loop is 




(3.4) 
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both of which are small compared to Galactic scales. 

For comparison, the characteristic scale of the loops formed at equipartition is 

l7nax,eq = Oluteq = ac/14kpC ( . _ ^7n4 ) " (^'S) 

\4.7 X lO^yrs / 

Galactic scales are ~ 1 — 100 kpc and microlensing can probe the full range of loops 
generated during the radiation era plus the small end of the loops generated during the 
matter era. All this assumes au order unity. Generally, the internal velocities associated 
with loops are relativistic. 

We want to answer two questions: What is the probability for lensing a single source 
at distance R by a distribution of loops at a given instant? How does the probability 
grow with time? 

Consider a small loop of size / at distance r. It lenses an angular area fli ^ {9E'r)l/r'^. 
The probability that a single background source at distance R is lensed is the ratio of 
the lensed angular area to the observed angular area in the direction of the source. We 
find 

and JF is the overdensity of loops in the Galaxy (J-" = 1 gives the lensing probability in 
a uniform Universe. We will estimate JF below.) Assuming the ordering Ig < lmax,eq < 
R < auttoday we find for a homogeneous loop distribution 

where H is the first moment of the loop distribution scaled to fiducial parameters. This 
approximate analytic result is based on two joined powerlaws for radiation and matter 
eras. In the appendix we show 

H{x, y) = (0.63 + 0.37^ + 0.04 log x) (3.8) 

teq \ 1.35 X lO^Vs 



Note that as T^G^ — > that H varies. In particular, x oc y oc l/F/jGyU and Pl oc 
^y h/Tr. Extending the loop distribution to smaller sizes {fi 0), therefore, does not 
overcome the effect of the decrease in Einstein radius. 

A numerical evaluation of Pl in the A-CDM model is presented in Figure |l] for the 
basic parameter space in Gfi — Tn we will consider. The string network has au = 0.3, 
ax, = 10~^, (3 = 1.6 with normalization set by agreement with simulations T = 43.6 
(see Appendix). Here and elsewhere the numerical results we present are based on the 
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Lensing Pi 
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Figure 1: Probability that a hne of sight to a source aX R = 100 kpc is microlensed by 
intervening string loops in a homogeneous Universe {J- = 1). 

fiducial A-CDM cosmology; the approximate analytic results illustrate the basic scalings. 

Now consider an experiment that stares at a given source and looks for the doubling 
in brightness on account of the passage of a loop along the hne of sight. One loop sweeps 
out an area per unit time ~ cZ/a/3 where the numerical factor crudely accounts for 
velocity orientation effects. The rate of change in the solid angle is d^i/dt ~ clj^r^. 
The lensing rate is 

which yields the lensing rate per source per year for a homogeneous loop distribution 

VlOOkpc; V ttoday ; ^ '^^ ^ ^ 

By contrast to the situation for Pl, when TrG^ we find Rl oc ^l/F^j/i. 

The numerical evaluation of Rl for the same A-CDM model as described above is 
shown in Figure |^. 

We have estimated the rate Rl assuming that the loop is moving relativistically. 
This is generally the case for the internal motions of the loop about its center of mass. 
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Lensing Ri per yr per object 
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Figure 2: Rate of microlensing along a line of sight to a source at -R = 100 kpc by intervening 
string loops in a homogeneous Universe {J- = 1). 

The velocity of the center of mass of a loop accreted to the Galaxy will be the halo 
velocity as we describe in more detail in the next section. New sources are lensed at 
a rate {vh/c)RL with repetition ~ c/vh- 

4. Cosmology of Small G/i Loops 

The center of mass velocity of string loops within the galaxy is a key parameter in 
characterizing the lensing rate for two reasons. First, if the velocity is less than the 
characteristic escape velocity from the Galaxy Vh then loops will accrete and be bound 
to the halo. Second, if the velocity of small loops is much less than c then a single source 
is lensed multiple times. 

The initial loop velocity is determined by the interactions in the string network. 
Intercommutations between string segments can generate relativistic center of mass mo- 
tions for the newly formed loops. Our treatment of the loop distribution is based on 
the assumption that all intercommutations occurs shortly after a horizon-crossing loop 
is chopped up. General scaling arguments suggest that the interaction rate diminishes 
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rapidly. Without interactions, the center of mass velocity suffers cosmic drag.^ 

Assume that the initial center of mass velocity is = 0.1c at time ti when the loop 



is born. At later time t the center of mass velocity is ^5 

ajU) 
a{t) 



X 



(4.1) 
(4.2) 



Figure | shows the fraction fsiow of the first moment of dN/dVdl attributed to slow 
moving objects. "Slow" means the loop's center of mass velocity today < 300 km s~^ 
in the A-CDM model with the fiducial parameters. The first moment is dominated by 



f. dN 

Fraction I 1 

dldV 



with v<300 km/s 




Figure 3: The fraction of the length- weighted loop density distribution with center of mass 
velocity < 300 km 

objects moving slowly enough to bind to the current halo. 

The scale of the gravitational potential induced by a low tension string is small 
compared to that of the amplitude of perturbations entering the horizon ~ 10~^. At 



^We assume that the recoil from emitted gravitational wave radiation is small. 
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equipartition ^lioop ~ 0.0016a/G';u/2 x 10 and the mass scale of the loop component 
which contributes most to at that epoch is small (~ O.SM© {j^yq^Y)- We first 
treat the low tension strings as "test particles" with respect to the baryons and CDM 
perturbations. As the perturbations begin to collapse after equipartition the slow moving 
loops fall into the gravitational potentials which will eventually give rise to galaxies. 
The observed overdensity of material in the Galaxy is 

J over K^'^J 
i 'A/ pcrit 

~ 1.5 X 10^ (4.4) 

where we have used = 0.3, Hq = 70 km s~^ Mpc~^, and a galactic mass of 2 x lO^^M© 
in a characteristic size of 20 kpc ||33|. The size scale has been selected in consideration of 



the location of stars of interest for microlensing. An estimate for the local enhancement 
of loops is 

f slow f over (^■^) 

i.e. roughly the overdensity of the Galaxy itself. 

The above estimate for JF is a lower bound because the actual dynamical evolution 
of the loops is not dissipationless. At equipartition, a loop begins to accrete matter. 
A novel aspect of the structure formation in this scenario is that suitably small loops 
will radiate and disappear but leave behind bound sub-galactic clumps of matter. Such 
objects are of considerable interest in their own right but here we focus on the mass 
density of the smallest surviving loops of mass Mg in today's Galaxy. Such a loop 
accretes a mass ~ (1 + Zeq)Mg with Zgg ~ 3570 (in the approximation of a flat, matter- 
dominated cosmology). This bound, composite object lives in the larger-scale growing 
galactic potential and it will be dragged by dynamical friction towards the center. As 
an example, consider the case for Gp = 2 x 10~^°, where Mg ~ 1.7 x lO^M© and clump 
mass ~ 6.4 x IO^Mq. Let the clump move on a circular orbit in an isothermal halo with 
rotation velocity 220 km s~^. The radial drift takes it from ~ 36 kpc to 20 kpc over a 
Hubble time [|^. The net motion with respect to the halo increases the density of such 



loops within 20 kpc by at least an additional factor ~ 1.8. The net radial drift for a loop 
of length / scales oc (/i/)^^^ and the drift for the smallest surviving loop scales oc jjy/Tji. 
These dissipational effects mean that loops, like baryons, may be over-concentrated with 
respect to the dark matter. The effect on the smallest loops is important if Gfi > 10^^'^. 
The distribution of loop size will also be tilted within the Galaxy compared to dN/ dVdl 
on account of the dissipative processes. 

We make a conservative estimate for T in our numerical calculations by ignoring 
these dissipative enhancements to the loop density. 

5. Practical Lensing 

There are a number of characteristic timescales relevant to experimental detection of 
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lensing. For the lensing itself, the characteristic duration of the event is 



A/ — 



c 



= 6.3 X lO-'sec — '—-r (5.1) 

VlOOkpcy V2 X 10-10/ ^ ^ 

The characteristic time for the (smallest) loop to oscillate is 



tosc ~ f ~ 135yrs j (5.2) 

and this governs the repetition timescale. 

The observational timescales are 6ti, the time for an individual observation, AT, the 
duration of the experiment, and SAT = AT/N^ep, the characteristic interval between 
observations where Nj-ep is the number of times a star is visited over the course of the 
experiment. 

To make this concrete, we will consider the Gaia mission which will monitor — 10^ 
stars for AT = 5 yrs with fluxes exceeding a mission-defined broad band magnitude 
G — 20. On average, the Broad-band photometer will observe each star N^ep — 80 times 
but the interval between observations is not fixed; it will vary from ~ 30 minutes to 
AT/Nrep- An individual observation is Sti ~ 3.3 seconds. Let 6AT ~ ATobs/Nrep be the 
characteristic interval between repeat observations. The mission observes ~ 10'^ stars at 
a time (note N^^NrepSti » AT); all astrometric and photometric results for individual 
stars are derived by detailed analysis of the joint observations. 

The accuracy of the photometry depends upon the brightness of the source. The 
limiting magnitude G = 20 corresponds to objects with a range of usual visual magni- 
tudes y ~ 20 — 25. Individual observations for sources with V — 20, 21, 22 and 23, for 
example, have relative flux accuracies 2.5%, 5%, 10% and 26% respectively and are more 
than adequate to see a factor of 2 change due to microlensing at the limiting magnitude. 

A star like the Sun (type G2V) will be visible to a limiting magnitude V = 20.2 or 
a distance approximately 12 kpc. To make estimates of the number of lensing events 
that Gaia is capable of observing we need to account for the distribution of stars and 
the "efficiencies" with which detections can be made. We assume here that the stars 
observed by Gaia are uniformly distributed in space and concentrate on the detection 
efficiency 

fdet f mag f size f time (^•"^) 

where fmag is the flux limit, fsize is the source size cutoff (so that the angular size 
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O < Qe), and fume accounts for sampling in time. We adopt a crude model 



time 



if 5tiens < 5ti 

if 6h < 6ti,^s < SAT 

1 if 5 AT < 5ti^ns < AT 
if AT<6t 



lens 

Any lensing event which is shorter than the length of a single observation or longer 
than the length of the entire experiment cannot be detected; all remaining events longer 
than the average interval between observations will be seen but only a fraction of events 
shorter than it will be. The assumption is that the observational intervals are uniformly 
spread between the shortest and longest periods. A more accurate description of the 
schedule of observations will readily improve the model for fume- To find the detectable 
lensing rate, we recalculate Rl with fdet to give i?^'**". The expected number of events 
is N^"'"" = Rl'^'^'N^AT. 

The expected number of lensings detectable by Gaia are shown in Figure ^ by 
red/blue contours. These include all efficiency factors (magnitude, size and time) with 
each contour labeled by the log of the number of detectable lensings. Note that the ex- 
pected number of detectable events decreases as Gfj, — >^ at fixed (for Gfi < 10~^), a 
consequence of decreased detection efficiency for short duration events {fume 0). The 
sharp cutoff as Gfi ~ 10~^^ is where the angular size of the star becomes comparable to 
the Einstein angle {fsize 0). The number of events depends on the path length (set 
by the flux limit fmag) and the intrinsic loop density. 

To illustrate the important effect of size and time cutoffs, we also display the expected 
number of lensing events if detections were limited only by flux considerations (green 
dashed contours). Note that the expected numbers increase as Gfi 0. These curves 
give some guidance on how alternative experiments might fare by adjusting Nrep, iV* and 
AT. 

An experiment with duration AT < tosc will see at most one lensing event. Con- 
versely, if AT > tosc the experiment has the capability of seeing multiple lensing events 
from the same source. If Gfi < 7 x 10^^^ multiple lensing may be observed by Gaia. In 
any case, a potential lens can be re-observed by other means even after the mission has 
ended. 

The figure has been constructed for a "typical" star within the Galaxy. The ob- 
jects included in Gala's catalog will depend on many factors - stellar type, interstellar 
obscuration, crowding, etc. - and these will alter the effective cutoffs and the number 
of events. For example, bright halo giants will be visible within approximately 41 kpc 
(type GOIII; radius ~ 3.3Rq). Compared to the more numerous G2V class, the lensing 
rate per star should be ~ 14 times larger on account of longer pathlengths and longer 
lensing durations. On the other hand, the cutoff at small Gn will be comparable to G2V 
because the angular sizes are comparable at the limiting distance. 
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Figure 4: Log of the number of microlensing events detectable by a Gaia-like mission looking 
at 10^ sources over 5 years. In this plot all sources are assumed to be solar mass stars (G2V) 
within 12 kpc. The string loop density is based on the fiducial powerlaw model evolved in a 
A-CDM cosmology. The local string overdensity is J'^ = f slow f over as described in the text. The 
red/blue lines include magnitude, time and angular efficiency factors; the dashed green lines 
include only the flux cutoff. 

At the flux limit, a star's angular size oc where Tg// is the star's effective 

temperature. Hotter stars can be seen further and have smaller angular size at the 
limiting distance. Consequently, the condition that the Einstein angle exceeds the stellar 
angular size implies a minimum tension detectable by an optical microlensing experiment 
oc The inclusion of stars hotter than the sun may provide sensitivity to tensions 

beyond the cutoff indicated in the figure. Of course, the number of such stars is an 
important consideration: hot white dwarfs, for example, would probe very low tensions 
but are not visible to large distances. Likewise, the duration of the lensing event should 
exceed the timescale of a single observation to measure the doubling in image brightness. 
At the flux limit, the minimum tension oc where L is the stellar luminosity. The 

inclusion of stars more luminous than the sun may provide sensitivity beyond the cutoff. 
An accurate assessment will depend on the make up of Gala's catalog, a subject to which 
we hope to return. 
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The Broad-band photometer measures 5 colors during an observation but this infor- 
mation is insufficient to allow identification and classification of stellar types throughout 
the HR diagram; it will be supplemented by photometry in 14 other bands from the 
Medium-band photometer, a distinct instrument on board. The mission expects to ob- 
serve ~ 18 X 10^ variable stars p2[ and the extensive color coverage is essential to 
facilitate identification and classification of known variables types (pulsating stars, erup- 
tive variables, etc.). This circumstance will prove a great benefit if string microlensing 
is sought. String microlensing is almost achromatic (the Kaiser Stebbins effect gives 
relative frequency shifts order ~ G^'~f{v/c); since Gfi is small the shift is small for mildly 
relativistic loop motion) whereas variable stars typically show color changes. The key 
point is that a cosmic string microlensing event should be an achromatic change in bright- 
ness by a factor of 2. The situation is qualitatively similar to that faced by ground-based 
microlensing experiments which anticipate an achromatic change in an a priori known 
form for the time-dependent amplification factor. 

Finally, we note that the Broad and Medium-band photometers have their own time 
sampling characteristics. Depending upon the detailed scheduling of scans, there may 
be an improved detection efficiency compared to our simple estimate. 



6. Remarks 

When string tension is low the density in loops exceeds that in long strings because 
the timescale for gravitational damping becomes long. The total string energy density 
today is dominated by the loops originally formed during the radiation era. This pre- 
ponderance of loops over long strings motivates consideration of the cosmology of loops. 
The evolution of these objects turns out to be completely distinct from that of the long 
strings originally studied with Gfi ^ 10~^. 

The loops damp by cosmic drag and are utterly inconsequential in terms of structure 
formation on scales much larger than the galaxy. Instead, they fall into the potential 
wells created by cold dark matter and baryons after equipartition. The Galactic halo 
number density of loops is enhanced over the universe's average value by at least the 
Galaxy's overdensity ~ 10^. Dissipational effects which depend upon the size of Gfi may 
further increase the overdensity. 

The loops within the Galaxy can be observed by microlensing using stars as sources. 
To the extent that stars are point-like and bright low tension strings can be observed. 
Simple estimates of a microlensing survey based on the capabilities of the Gaia mission 
suggest that many such events may be detected. If loops are moving at halo velocities 
the lensing of a given source should repeat ~ 10^ times. 

Detailed observations of such a lensing source will have much to tell us about string 
tension and the number density of loops. 

We recognize many gaps in this general story 



- 18 - 



• What processes are responsible for producing the T — 1 (average, unclustered) 
distribution of loop sizes? When does the rate of interactions between loops and/or 
long strings become negligible? When does cosmic drag predominate over velocities 
induced by intercommutation? 

• How do loop-loop interactions develop as the galaxy grows? Does intercommuta- 
tion end up ejecting loops from the galaxy? How does the distribution function of 
loop sizes change on account of the intcrcommutations in the dense environment? 
How small arc the smallest loops? Is the gravitational lifetime significantly al- 
tered? Can one use this to measure or learn about p7 The loop size distribution is 
potentially measurable via microlensing and locally generated gravitational waves. 

• What are the astrophysical constraints on massive loops moving around galaxy? 
How are loops distributed within the Galaxy? 

• What string and loop parameters can we deduce by observing repetitive lensing? 
Does small scale structure on the loops increase the microlensing repetition rate? 

• Incorporate lensing over radii ranging from larger than the loop size to smaller. Is 
the Kaiser-Stebbins effect detectable from a known lens? 

• An accurate calculation of Gala's capabilities for string microlensing requires in 
depth analysis of the time-based observing strategy for both the broad and medium 
band photometers. 

• Can one design an experiment tailored to doing a better job for detecting microlens- 
ing than a mission like Gaia? It is worth noting that the Gaia sample includes pho- 
tometry much more accurate than what is needed to detect microlensing. Could 
one generate a larger catalog with more poorly determined flux measurements still 
adequate to detect string microlensing? Would there be enough information to 
rule out variable stars? 

• How can one optimize the choice of AT, Nrep and 6ti to investigate particular 
ranges of G/x? 

Detecting cosmic superstrings and measuring their tension can reveal fundamental 
information about superstring theory and elucidate the large-scale contents as well as 
the remote evolutionary history of our own universe. An observable, nearby source of 
strings to study promises to advance, qualitatively and quantitatively, these goals. 
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A. Appendix 

We start with a form generally consistent with the assumption of a scaling network 

r^f{x) with X = l/t (A.l) 



dN{l,t) _^^5 



dldtdV 

where / and V are measured in physical (not comoving) lengths. We make the assumption 
that the chopping of long strings all the way down to the smallest loops occurs in a time 
short compared to the expansion timescale. Of course, on the Hubble scale this won't 
be a very good approximation but it should get better at smaller scales and thats what 
is of most interest here. Loop length I is the time-averaged physical length associated 
with the loop when the loop itself is small compared to the scale of the horizon. A 
cosmological simulation is required to establish the form for / which describes loops 
with sizes comparable to the scale of the horizon. 



Motivated by the numerical simulation in Ref. l23, 28], we take the loop production 



to be given by a power law distribution with upper and lower cutoffs 

/(x) = Ax~'^ for ai < X < au (A.2) 

and zero otherwise. In addition to the upper limit, we impose a lower limit on the 
the range of loops. This is a manifestation of the cutting up process and not due to 
gravitational wave emission. In this approximate description, dynamical processes chop 
up strings to give the loop spectrum instantly. Gravitational radiation acts once the 
spectrum has been created. 

The scaling network is characterized by some inter-string distance d{t) = 'jst, defined 
so that the density in long strings is poo = p/c?^- Conservation of energy then gives 

dldVdt dt ^ ' 

or 

xf{x)dx = 2^ (1 - {v^)) (A.4) 

where (f ^) is the square of string velocity averaged along the length of long strings. 
The energy balance between long strings and the loop distribution gives 

A= ^(V<"'»P-g) ^ (A.5) 

2 2-/3 /-I / / \2-/3^ ^ ^ 
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Numerical simulations derived A = 82±2, ai = 10 ^, au = 10 ^ and P = 1.63±0.03. 
From these we deduce the numerical value for 



T = 




43.6 



(A.6) 



which we will hold fixed even as we consider models with different characteristic values 
of G/i, aL, au, P and different expansion dynamics n ^ 2/3). 

It is generally thought that a loop of length / decays by gravitational radiation in 
time Z/(rG/i), where F is the dimensionless decay rate of order 50. We write our model 
with two generalizations of the canonical description. 

• F/j: We introduce a distinct damping rate for loops, F^?. It is possible that the 
average loop damping rate differs substantially from F because small loops, large 
loops and horizon-crossing strings intercommutate. 

• r: Lifetime in the network might be a nonlinear function of /. This could arise 
because of the nontrivial interaction of intercommutation and damping. 

If lifetime is linear in length then loops of length / born at time t' and observed at 
time t must have 



Our model assumes that loops are unaffected by gravitational wave damping until they 
reach the end of their life at which point they are abruptly removed from the population. 

The loop density at time due to loops produced during the interval X = (to, ^i) is 
found by integrating the production rate density production of loops: 



for t > t' 



(A.7) 



dVdl 



dN 




dVdldt' 



dN 



v{t')e{t2-t' <t{i)) 



(A.8) 



Here V{t) is physical volume. Inserting the production rate density 




(A.9) 



where 




(A.IO) 




(A.ll) 
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A.l Power law / and a(t) 

During I assume that the scale factor a oc t", the volume V oc t^^, and f{x) = Ax~^ for 
ai < X < au (and otherwise). The loop density at t2 is 



dN , ^ A 



" J.* "I C 

f* 



(A.12) 



where 

( = 4-3n-/5. (A.13) 

For typical values of (3 and n we have C > 0- Now the density at t2 depends upon how 
V varies beyond X. If n is constant over the entire interval to to ^2 then Q = 1. If loops 
are born in the radiation-dominated epoch [n = rii = 1/2), which extends from t ~ to 
teq ~ 4.7 X 10'^ yrs, and are observed in the matter-dominated and/or lambda-dominated 
epochs we have a non- trivial Q. For the ordering h < < t2, G = V+{teq)t^''^ /V+{t)tl'^^ 
where V+ is the volume for t > t^q. If a{t) oc t"^ for t > t^q then g = (t/teg)^^"'""'^- This 
reduces to ^ = (t/teq)'^^"^ for rii = 1/2 and n2 = 2/3. 

It will prove useful to rewrite the general result in the following form 

g^f^{^,] (—,) (i-(^) 1 (A.14) 



dVdl 




(A.15) 



Figure || presents an estimate of the loop density today obtained by adding the loops 
created during the radiation and matter-dominated eras. Specifically, we take n = 1/2, 
to = 0, ti = teq followed by n = 2/3, to = teq, ti = t2 = ttoday For lack of more precise 
numbers we take ai = 10""^, au = 0.3, T^Gfi = ac = 10~® (i.e. lifetime linear in /: 
r(/) = l/ac), P = 1-6, T = 43.6 during both epochs. These choices imply A/i = 40 
and A/2 = 91. The background cosmological model has te = 4.7 x 10^ yrs and the 
current epoch ttoday = 1-35 x 10^° yrs. The estimate of the number density has a cutoff 
at / ~ 10~^ from gravitational damping; the very slight kink at loop size / ^ 10^^ is 
roughly equal to auteq i.e. of order the horizon at equipartition. 

Two straight lines are included in the figure and are evidently close approximations 
to the slope of dN/dVdl over an intermediate regime. The slope is determined by 
the expansion rate. Assume a pure powerlaw expansion with to = and ti = ttoday, 
ai « au and C > (e.g. f3 = 1.6 and 1/2 < n < 2/3 gives 0.9 > C > 0.4). 
Define the lengths /^m = c^Lt, Imax = ctc/^ and Ig = act. For I large compared to both 
the gravitational damping and the minimum chopping scales ^o ^ ^ ^ small 

compared to the horizon / -C Imax we have 



dVdl /4 V t 
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Figure 5: The number density of today's string distribution based on two powerlaw semi- 
analytic estimate (long dashs) and ACDM model (solid line). Two straight lines (short dashs) 
with slope -2.5 (radiation) and -2 (matter) are included for comparison. 



and 



^=?<i^ (A.17) 



This form shows that the slope of the loop distribution at intermediate scales is com- 
pletely determined by the rate of expansion: dN/dVdl oc /"^'^ for n = 1/2 and oc for 
n = 2/3 as previously demonstrated by Vanchurin et al. 



Figure ^compares the distributions at t = ttoday and at IQ'Htoday The gravitational 
wave damping destroys the small scale loops while the expanding horizon continuously 
gives rise to larger scale loops. The expansion decreases the number density. The inter- 
mediate form, supplemented by the appropriate will generally be sufficient to describe 
the loop distribution for loops that cannot have decayed (/ > lg^2 where lg 2 = ^^(^2)) and 
smaller than the largest loop formed (/ < lmax,i where lmax,i = c^uti), presuming the 
ordering lg^2 < I < lmax,i- Of course, the full expression is needed to describe the exact 
form near small and large / and to deal with situations in which the intermediate regime 
does not have enough time to form. 

We are primarily interested in understanding whether various integrals are domi- 
nated by the small or large scale part of the distribution. Consider the density distribu- 
tion of Z*^; (5 = gives number density, Q = 1 gives length (needed for energy density 
and lensing estimates), and so forth. Using the intermediate form for lg^2 < I < lmax,i 

Clearly, Q < 3 — Sra is dominated by small scale loops (Q < 3/2 for n = 1/2; Q < 1 for 
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Figure 6: Number density of tiie string length distribution at two epochs ttoday and 10 '^Uoday 
in the two powerlaw model. 

n = 2/3). Assuming auti > act the integrations give the relatively simple forms 



dN TV 2 L/^ - ./f if n=l/2 

l^it-M)dl ^ { V *^y (A.19) 

dVdl I log (^^^ if n=2/3 

The total integral is formed from the sum of loops created during the radiation 
and matter-dominated epochs {Qi = {teq/tY^"^, Q2 = !)• We assume (for lack of any 
more complete information) [3i = P2, Ciu,i = «{/,2 and aG,i = aG,2- The final result for 

t ttoday is 



l^dl ~ i f 2A/-, - lU Af, log f ^ 1 1 , (A.20) 



dVdl \ \V ctct J \0!G 

If acttoday C(uteg the loops Created before t^q dominate and 



dN „ 2Afl [out eg 



The total energy density is 



l^TT^dl^^J^^^ (A.21) 
dVdl \ aat ^ ^ 



Figure ^ plots t^dN/ dVdl so that contributions in different logarithmic intervals of / can 
be directly compared. Evidently, for the parameters adopted the small loops formed 
during at t < tg dominate pioop- This result (and other qualitative results alluded to 
above) will be altered if ac/au > teq/Uoday ~ 3.5 x 10"^. 

For a fiat, matter-dominated cosmology (i.e. ignoring A) we find that 

ni^op = QTrGf, (^2^, -l^+^2 log ) . (A.23) 
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Figure 7: Today's string number distribution weighted by length; Hues as in previous plot; 
straight (short dashed) lines have slope -.5 (radiation) and (matter). 



A. 2 Analytic Estimate 

Define two dimensionless parameters, scale to "typical" parameters, as follows: 



VO.3/ 

Xo = 3 X 10^ 



X 



t 



eq 



XqJ V4.7 X 10 



1.35 X lO^Vs 



t 



today 



Yo = 104.4 



which gives 



dVdl 



2(2 -/?)T 



U2 17.22 + log 



X 



au 



2-/3 



(A.24) 
(A.25) 
(A.26) 
(A.27) 



( = 4 - 3r2 - /3 

Adopting numerical parameters approximately consistent with the VS simulations 



(A.28) 
(A.29) 

(A.30) 



au = 0.3 
tti = 10"^ 
f3 = 1.6 
T = 43.6 



(A.31) 
(A.32) 
(A.33) 
(A.34) 
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we find 



and 



/ 



Ni = 40.4 (A.35) 
M = 90.9 (A.36) 



nioops - 8.5 X 10-^ (2^^^) H{X, Y) (A.38) 



H{X, Y) - I 0.63 + 0.37^^ + 0.04 log ) ) (A.39) 

where X and Y depend upon au/TnG/j, and cosmological timescales. 
A.3 A-CDM 

For a more realistic cosmology, consider a A-CDM model. Let ^2^,0 = 8.4 x 10~^ ^m,o 
0.3 and Qa.o = 1 ~ ^r,o ~ ^m,o- The relationship between a and t is given by 



/° 

^0 



, da _ _ = r (A.40) 

where r(a) = Hot{a). We choose a = 1 today {t{l) = ttoday, -f^o = r{l)/ttoday) or 
t(a) = UodayTifij/Ti^) s-nd l^(a) (X a^. Assuming that all parameterized quantities 
au, OiL, 0(G, 7s) s-re constant, we evaluate the loop density at Uoday numerically 



dN A 



r dt't''^-^V{t') (A.41) 



dVdl If'Vittoday) Jt*, 

where 

tl = min (^ttoday, (A.42) 

^0 = max (0, ttoday - — , — ) • (A.43) 

These results are included on the plots; they are the solid, smooth lines. They agree 
fairly well with the simple, two-powerlaw approximation. 

A. 4 A-CDM: Numerical Results for varying Gji, au, (Xl and P 

We choose a fiducial set of parameters: /3 = 1.6, TrG/j, = 10~^, au — 0.3 and az — 10~^. 
We then varied a single parameter over a wide range to explore how the loop distribution 
function changed. 
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Figure 8: Loop distribution for TjiGfi = 10 to 10 ^ in steps of 10^ for standard case. 
Slight offsets to expose separate curves. Lower TjiGfi extend to smaller loop size. 

Effect of variation of Gfi on the loop distribution is illustrated in Figure Decreasing 
the tension extends the distribution to smaller size scales because gravitational radiation 
is less rapid. 

Effect of variation of (3 on the loop distribution is illustrated in Figure ^ Smaller 
f3 concentrates loops at the largest scale relative to the horizon. These large loops live 
longer with respect to gravitational wave damping modestly increasing the abundance 
of the smallest loops. 




Figure 9: Loop distribution for /3 = 0.6 (red), 1.6 (green) and 2.6 (blue) 



Effect of variation of au on the loop distribution is illustrated in Figure |I0| Smaller 
au limits the size of the loops at all times and leads to more rapid decay. 

Effect of variation of a/, on the loop distribution is illustrated in Figure |n|. Larger 
ai increases the size of the loops at all times and leads to less rapid decay. 
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Figure 10: Loop distribution for au = 3 x 10 ^ (red), 3 x 10 ^ (green), 3 x 10 ^ (blue), 0.3 
(black) . 
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